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Abstract 

The standard way of evaluating residues and some real integrals through the residue theorem 
(Cauchy's theorem) is well-known and widely applied in many branches of Physics. Herein we 
present an alternative technique based on the negative dimensional integration method (NDIM) 
originally developed to handle Feynman integrals. The advantage of this new technique is that 
we need only to apply Gaussian integration and solve systems of linear algebraic equations, with 
no need to determine the poles themselves or their residues, as well as obtaining a whole class of 
results for differing orders of poles simultaneously. 
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I. INTRODUCTION 



In a textbook on complex variables we may find real definite integrals of the type 

P dx f dec 

' ^ Jo Jo (x2 + l)(x2+4) ^ ^ 

to be evaluated using the contour integration in the complex plane, making use of the 
Cauchy's residue theorem. This is, of course, a simple exercise in complex analysis, and 
residue theorem is a powerful tool to handle such integrals. However, if we were actually to 
evaluate them we would do it separately, each integral in its turn, with its own residues of 
poles summed over to get the final answer. Not so if we use the NDIM technique, as we shall 
shortly see. In NDIM we can integrate both integrals at the same time, and more, without 
having to do it one by one separately; we simply need to evaluate the general integral 

T(m,n,p) ^ dx (X^)"^ 

and then work out the result for a particular set of numbers m,n,p = 0, 1, 2, 3, of which 
Ji and I2 are specific examples. 

The concept of negative dimensional integration can be best seen in the following D- 
dimensional Gaussian integration: 

G^/Ae--^(^° (3) 

On the other hand, expanding in power series the exponential function in the integrand 
of (0) above, we have 

G = I d^xe 

^ n 



— [X 

nl 

n=0 



d^xJ2H 

n=0 

°° n r 

Now, comparing (jS)) and we conclude that for the equality to hold we must have 

X ^ 1 d''x{xy = {-l)-n\ iV^f 5^^o^o (5) 
and since, by construction n > 0, one is led to assume that D is a negative valued dimension. 



This non-trivial result for a D-dimensional integral with pure quadratic integrand ele- 

n 

vated to a given power was first pointed out by Ricotta and Halliday in 1], which differed 
somewhat from all the previous considerations, starting from dimensional regularization 
scheme developed by 't Hooft and Veltman [2], where all such integrals were set to zero 
straighforwardly. Since Feynman loop diagrams lead to Feynman integrals in D-dimensions 
in the dimensional regularization scheme, integrals evaluated with the help of (0) in negative 
dimensions, must be brought back by analytic continuation to the realm of positive dimen- 
sions D. Since the concept of negative dimensional integral arises as a consequence of the 
positiveness of the polynomial powers in the integrands, the mentioned analytic continuation 
is achieved by employing a property of the Pocchammer's symbols bearingthese parameters 



labeling powers of the polynomial expressions in the integrands, namely 

(a)_i. 



i 



(6) 



(1 - a)k 
where 

Of course, here we strict ourselves to one-dimensional integrals such as (j2)), so that (jSj) 
becomes 

Id=i = j dx{xy = {-iy^n\ ^ 5^^i_^^ (7) 

and all the NDIM formulation developed for D-dimensional Feynman integrals can be applied 
to evaluate one-dimensional integrals of the type Ii and l2- 



II. REVIEW OF CAUCHY'S RESIDUE THEOREM APPLICATION 

Consider, for practical example, the integral 

dx _ 1 /'+°° dx 

The second integral on the RHS of the above represents an integration along the real axis 
of the function 

^(^) = = ^z + i){z-i) 

which is a function of complex variables with simple poles at z = ±i. 

Let Cr be the semi-circle with radius \z\ = R with i? > 1 as shown in Fig. 1 below 
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Im(x) 




By Cauchy's theorem, we have 



R 



dx f{x) + / dz f{z) = 2Tci Ki 

R JCr 



(10) 



where ki is the residue of f{z) at the (simple) pole z = which, of course, can be calculated 
easily by 



Ki = {z-i)f{z)\, 
Therefore, for i? > 1, we have from ()10|) 



1 



z + i 



1 

2i 



R 



dx 



Now, \z\ = R when z is on the semi-circle C/j, so that 

\z^ + 11 > UP -1 = -I 



TV 



dz 



z^ + 1 



and therefore 



dz 



< 



\dz\ nR 



for i? — > oo 



Then, taking the limit. 



R 



dx 



or 



lim , 

R^oo J_^ X'^ + 1 



r+co 



dx 



x^ + 1 



TV 



dx TT 



x^ + 1 2 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 
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III. NEGATIVE DIMENSION INTEGRATION METHOD 



On the other hand, let us take the Gaussian generating funtional of the negative dimen- 
sional integration, namely, 

r+oo 

1= / (ixe-"("'+i) (17) 



where we note that the argument in the exponential is the denominator of the integrand in 
(jH)) multiplied by a real, positive parameter a which is chosen as a converging factor for the 
integral. Of course, this is the standard Gaussian integral, whose result is 

oo , X _1 

If we project out in power series the integrand of (fT7|) before the integration is done, we 
have 

X=$^(-r^/NDiM(n) (19) 

n=0 

where the negative dimensional integral /ndimI?^) is given by 

r+oo 

hDiM{n)= / rfx(2:2 + l)- (20) 



Note that for n = —1 this is exactly the integral (jH)) we want to evaluate. Observe, 
however, that the negative dimensional integral is more general than the one we have in (jH} in 
the sense that the exponent n is not fixed. Note however that (j2Uj) is only defined for positive 
n, so that in order to get the result for (jHI) we need to make an analytic continuation to 
negative values of n, a process whereby the integral gets defined into positive dimensionality 
in (1201). 

Comparing the two series expansion (jlSp with ()19p . we observe that in order to both 
series be equivalent, we must have n = j — |, so that 



,1 r(n + l 

-TT) 2 



(n + 1) 



r(n + i + l) 

(21) 



2 

'-7r)2 
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Now, analytic continuing (AC) the exponent n into negative values using the property 
I, we have 

.(AC) , , f^°° dx 1 , , 1 r(-n - i) 

Now, substituting n = — 1, — 2, — 3, we get 

j(AC) , 2^ _ _ .24) 

^NDImI - (^2 + 1)2-2 ^^^^ 

j(AC) , 3X _ _ 3vr , . 



and so on and so forth. Of course, for the original integrals, which range from [0, oo), the 
corresponding values are half of the values quoted above. 

A more interesting case is the evaluation of the other integral, namely, 

' = X (x2 + l)(x2 + 4) = 2 (x2 + l)(x2 + 4) ( ^ 

Evaluation of this integral by the residue technique is similar to the previous one where 
now we have two poles in the upper hemisphere, z = i and z = 2i, so that the contour Cr 
is such that its radius must be \z\ = R, R > 2, before taking the limit R oo. The result, 
after summing the residues of both poles is 

= -2in< + — I = - <^ + _ k = (27) 

'2 \6z3zj2\33j6 ^ ^ 

We have written down the explicit contributions of each residue, the first one correspond- 
ing to the residue at the simple pole z = i and the second one to the residue at 2 = 2z, to 
show that these residues correlate with each of the "basis" (linearly independent) solutions 
(with a word borrowed from the language of basis vectors in a vector space) in NDIM. 
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The Gaussian generating functional of the negative dimensional integral to this case is 

p+oo r+oo 



— oo 

1 

712 



e-/5-47 



(« + /? + 7)5 



1 00 or 1 

(a + /3 + 7)^^^o 

a+b+c=-^ , , 

1 . , 2 aob+r c+s 

= ^^r(-) y (-r-4- "f„ ,1, (28) 

^2^ ^ ^ ^ a!6!c!r!s! ^ ^ 

where in the last line of the above we have employed the standard multinomial expansion 
for (a + P + 7)~5. 

On the other hand, direct expansion in power series of the integrand yields 
J = y (-)^+'+-^^ / dx {x^ix' + l)\x' + 4)™ 

= E i-y^'''"'yjrT^NBmU,l,rn) (29) 

j,l,m=0 ■' 



Comparing ()28|) with ()29|) we see that we must have 

j = a 
/ = 6 + r 
m = c + s 



a + 6 + c = (30) 



Also, when we compare the two series, the one in (jzHj) has five summation indices with 
one constraint, whereas the other one ()29p has three summation indices. Then our result 
for /ndim() i) ^5 '"^) will be given in terms of 5 — 3 — 1 = 1, that is, a single summation 
index. However, since we have five indices in one and three in the other, with one constraint 
equation among them, one can in principle have the following combinatorial possibilities: 

5 



5C^4 = ^ = 5 (31) 



for the remaining series index. 
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Letting r be the summation index in the result, we have then the following conditions: 

a = j 
b ^ l-r 



c = -j -l-- + r 
s = j + l + m+ - — r 
so that the negative dimensional integral will be given by 



X 



r=0 



j! (/ - r)! (-j - / - i + r)! r! (j + Z + m + I - r)! 



r(i)r(i + m) 



r(i-j-/-i)r(i + j + / + m + |) 



X 



E 

r=0 



(1 + 0-r(l -j-l- |)r(l + J + / + m + |)-r4^r! 



(-0. {-j - 1 



2 r=0 



-TT 2 



(i)_,_;(l + m) 




(32) 



(33) 



Analogous calculations can be performed for s, 6, c indices, whereas for the index a there 
is no solution except the trivial one (This is the case when the determinant of the system 
of linear equations vanishes). For the remaining non- vanishing determinants, we have then 
the solutions: For s remaining summation index, we have the conditions 



and the solution yields: 



'NDIM 



a 
b 
c 
r 



= J 



-j -m - - + s 
—m — s 



j + Z + m+ -- s 



-tt) 2 



X)-3-m (1 + Oi+TO+i 

(-7r)5 



{\)-3-m (1 + Oj+m+i 



E 

s=C 
2F1 



{-m)s {-j - I 



m 



(1 - J - m- i). 



-m 



-j -l-m-\ 



1 - j - m - I 



s\ 



(34) 



(35) 
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For b summation, we have 



and the solution yields: 



rib) 



'NDIM 



{j,l,m) = 



a = J 
r ^ l-b 



-TT 2 



j + m + - + b 
•' 2 



5 4i+'i 



Finally, for c summation we have 




(36) 



(37) 



a 



s — m 



b — —j — t: — c 



and the solution yields: 



r(c) 
'NDIM 



(i)-.(i + 



-m 



(38) 



(39) 



Now, observe that we have four non- vanishing solutions arising from the solving of systems 
of linear equations, which here clearly comes in pairs when we look at the argument of the 
hypergeometric functions in the results. The solution for the integral is then the sum of the 
pairs with the same argument (linear combination), namely. 



InmM{j,l-,m) = A2Fi{a,b;c\z ^) + B 2Fi{d,e; f\z ^) 
= A' 2F,{a', b'; c'\z) + B' ^F.id', e'; f'\z) 



(40) 
(41) 
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where the coefficients A, B, A' and B' are given by 

A 



B 
A' 
B' 



(1 + m)^.+;+ 
(-7r)5 4™ 



(i).,(l + m)^.+ i- 



and the hypergeometric function parameters and variables are: 



(42) 



a 


b 


c 


d 


e 


f 




-I 


—j — I — m — ^ 




—m 






1 
4 


a' 


y 


d 


d' 


e' 


f 


z 


—m 


-i -l-m-\ 


I- ] -m-\ 


-I 


3 + \ 


1+j + m + l 
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Before proceeding, let us demonstrate that the two sets of solutions, namely the primed 
and unprimed ones are totally equivalent. To do this, we employ the following analytic 
continuation property of hypergeometric functions jsl 



2Fi{a,b;c\z) 



(1 - b)a (c)_„ 



2Fi{a,l + a-c;l + a-b\z'') 



+ 



2Fi{b,l + b - c;l + b - a\z 



1^1 - a)b [cj^b 

to one of the basis solutions, say of the unprimed set 



\arg{-z)\ < n (43) 



iFi{-l, -J -l-m 



2' 



J -I- 



- I 

2 I 4^ 



to get the primed result. Of course, we could have used the transformation property above 
to the other basis solution 2Fi{—m, i + + / + | |^), and we would get the same primed 

result. We need to apply only to one of the basis solutions, since the transformation property 
above referred to cannot produce neither new nor any more than two linearly independent 
hypergeometric functions. 

Now, we need to analytic continue the results to negative values of exponents and positive 
dimension. Observe that the result for /ndim contains two factors: One is the coefficients, 
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given by ratios of gamma functions and the other is the functional part, given by the hy- 
pergeometric funtions. For the coefficients, which contain ratios of gamma functions given 
in terms of Pocchhammers symbols, we employ (jH)), and for the functional part, just let 
the exponents go to negative valued parameters. Note that we need to be aware of which 
exponent should be continued to negative values Then, our final result for the integral 
reads: 

^muUJ^rn) = A^"" ,F,{a,b;c\z-') + B'''' 2Fi{d,e-J\z-') (44) 
= A'^"" 2F,{a', b'; c'\z) + 5'^^ ^^^(^/^ g/. (45) 

where 

5AC ^ ^l4-(l),H).^_i; 

B'AC ^ ^\ 4i+m+i ^^^^ (-m)_._i. (46) 

One interesting thing about the NDIM techmology is that it allows us to write the correct 
answer in as many equivalent ways as it is possible to do. For the case in question, we have 
two equivalent answers, namely, the unprimed and primed answers. 

So, for particular values of the exponents, say, i = 1 and I = m = —1, which is the case 
for the I2 integral mentioned in the introduction, we have: 

^NgiM(l, -1, -1) = {4-^ 2F,{1, i; III) + 4-1 (l)i (1)_3 I; ||1)} (47) 

= vr^ {(l)-i 2^1(1, |; ||4) + 4i (|)i (1)_3 ,Fi(l, |; ||4)} . (48) 

Finally, using the fact that 

2F,{a,b;b\z) = {l-zr- (49) 
the two results (jTTj) and (jlH)) coalesce into one: 



dM(l,-l,-l) = -| + ^ = ^ (50) 



so that each of the basis solution corresponds exactly to the residue of the poles at z = i 
and z = 2i. Finally, 

l2 = ll^Ema,-h-l) = l- (51) 
11 



Other particular cases, such as j = m = 0, / = —1, j — 0, I — m — —1 and j — 0, I — 
— 1, m — —2 can be calculated from the general solution, yielding. 



-^1 ~ 2 -^NDIm(O) 1)0)— 2' 

.(0,1,1) _ 1 7-AC ('n _i _n 

-'2 — 2 -'NDIMV'-') ^) 



2 -'JNUIMV"' -"-/ 12 ' 

Stt 
288" 



7-(0,l,2) _ 1 tAC (n _i _r,\ _ Stt 
-'2 — 2 -'NDIMI*-') -■-) ~ 288- 



Note that for all the cases where either / = or m = the general solution is such 
that one of the terms in the overall result vanishes because we have a term proportional to 
= 0, and the answer is then given by a single hypergeometric function. 

IV. CONCLUSIONS 

Using the NDIM technique, we evaluated some sample real definite integrals which may 
be calculated by the Cauchy residue theorem in the complex plane. The alternative method- 
ology here presented gives us the bonus in that all the generic exponents of integrands can 
be calculated at once, from where particular solutions can be drawn. There is no difficulty in 
the performing of the integration since the integration involved is of the Gaussian type and 
the technique requires only series comparison term by term and the solving of systems of 
algebraic linear equations resulting from such a comparison. We showed that the results for 
different variables are obtained simultaneously and they are equivalent to each other. More- 
over, for each set of basis solutions correspond the residue of a given pole. The strength of 
this new technique also can be envisaged in that simple, double or higher order poles can 
be evaluated all at once. 
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